Comment on "Saturation of the all-optical Kerr effect 



In a recent Letter [T], Bree, Dcmircan, and Stein- 
meyer calculated higher-order Kerr coefficients in no- 
ble gases using a generalized Kramers-Kronig (KK) re- 
lation applied to high order nonlinear processes. They 
found that the nonlinear index An for argon obtained 
from An = J2 saturates and goes negative in the 

intensity range / ^40-50 TW/cm^, well below the thresh- 
old for ionization. 

While a recent experiment has shown conclusively that 
there is no nonlinear index saturation below the ioniza- 
tion threshold [2], here we show on theoretical grounds 
that the calculation leading Bree et al. to their conclu- 
sion is incorrect. It errs by using multiphoton ionization 
(MPI) rates in a regime where perturbation theory breaks 
down. 

The KK relation used by [1] is 

nM^^v r{n + c.)^;^dn, (1) 

where is the coefficient for /c*'' order MPI with rate 
Wk = cFkl'^^ ^ ~ {VL + kuj)/{k + 1), and V denotes princi- 
pal value. The coefficients at used by Bree et al. (Eq. (1) 
of [T]) are obtained from the 7 — > oo or MPI limit of a 
recently derived ionization rate R [3], where, in atomic 
units, 7 = (2IpY/'^u} / Eq is the Keldysh parameter and 
where Ip is the ionization potential, uj is the laser fre- 
quency and Eq is the peak laser field. The rate R (Eq. 
(6) of ^) handles both the 7 > 1 MPI limit and the 
7^1 tunneling limit, and provides a reasonably accu- 
rate interpolation for 7 lying in between. 

In Fig. 1 we plot R for argon {Ip = 15.76 eV) for 
TiLu = 1.55 eV (A = 800 nm) as a function of intensity 
(red dashes). At low intensities R cx where M {—11 
for Ar) is the minimum number of photons for ionization, 
and at higher intensities R shows the characteristic roll 
off as the tunneling regime is approached. We also plot 
the total MPI rate w = X]fc=ii ""^fc (blue dots), where the 
contributions for fc > 11 correspond to above threshold 
ionization. Above ~ 70 TW/cm^ , w diverges badly. 
However, even as low as 5 TW/cm^, R and w differ by 
a factor of ~ 2. Thus we immediately see that the terms 
containing ak contributing to w are quite inadequate to 
describe ionization in the > 40 TW/cm^ region where 
the novel behavior of An is claimed. Thus, these are 
problematic for use in Eq. (1). Forcing a mathematical fit 
of w to i?, using the ak as free parameters, would require 
some of the ak to be negative for large k, an unphysical 
result. 

We now show how Bree et al. obtain saturation and 
negative excursion of An. In [1] it is stated that n2k > 
for fc < 9 and n2k < for fc > 10. It is the fc > 10 con- 
tributions which cause An to saturate and go negative. 
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FIG. 1: (color online) Ionization rates w (blue dots), Wk (solid 
black), and R (red dashes) calculated for Ar using the models 
in [1] and [3]. Here Wk is the fc**^ order MPI rate applicable 
for 7^1 (shown for k = 11), w is the sum over Wk from 
fc = 11 to 50, and R is the rate from [3] applicable from the 
MPI though tunneling regimes (7 ^ 1). 



First, we note that it is the < fl < uj portion of the 
KK integral (1) which can contribute negatively to n2k- 
In that interval, ^ in ak+i{^) ranges from u)k/{k + 1) to 
w. However, (7k+i{£,) = for ^ < ^thresh = Ip/[{k + 1)K\, 
as appropriate for a MPI coefficient [3]. Suppose k = 1. 
Then the KK integral would contribute negatively for 
uj/2 < ^ < ut, except for the fact that (72(0 — in this 
range. Thus n2 > 0. For the same reasons, n2k > for 
fc < 9. However, for fc = 10, the KK integral contributes 
negatively in the interval (10/11)ll! < ^ < cj because 
^ ~ ^thresh and o'ii{£,) > 0. Similar negative contribu- 
tions occur for fc > 10, where ^ > ^thresh- Thus, all 
negative Kerr coefficients n2k for fc > 10 are computed 
using cr/c+i(0 foi' £. ^ a. region for which 7 ~ 1.8 at 40 
TW/cm^ in Ar. In using Eq. (1) to calculate the terms 
n2k in An = '^n2jl-' , 7 has therefore been implicitly 
sampled from 7 = 10 (fc = 1) to 7 = 1.8 (fc = 10), a range 
over which an underlying theory of ionization should en- 
compass the transition from MPI to tunneling. Instead, 
Bree et al. have used MPI exclusively, as appropriate 
for a perturbative expansion. It is not a coincidence that 
n2k < for fc > 10: this is the A = 800 nm MPI threshold 
for Ar, where crii(cLi) > 0. Essentially, the saturated and 
negative intensity dependence of An obtained by Bree et 
al. is an artifact of a perturbation expansion in a laser 
intensity region where it is far from appropriate. We note 
that Bree et al. justify their use of rates in the MPI limit 
by quoting 7 = 1.62 at 50 TW/cm^. While this might be 
an adequate assumption if the rates are to be used only 
for order-of-magnitude estimates of ionization yields, it 
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is clearly incorrect for calculating nonlinear coefficients 
using KK theory. 

Finally, we note the remarkable coincidence that an 
experiment [1] and a theory jl] both show apparent sat- 
uration of the Kerr effect in the same intensity range 
for completely different reasons. In the experiment the 
saturating and negative phase shift was likely caused by 
diffraction from a plasma grating [5] ; here it is caused by 
the extension of a perturbative theory of ionization to an 
intensity range where it breaks down. 
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